Indirect coupling between spins in semiconductor quantum dots 
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The optically induced indirect exchange interaction between spins in two quantum dots is inves- 
tigated theoretically. We present a microscopic formulation of the interaction between the localized 
spin and the itinerant carriers including the effects of correlation, using a set of canonical transfor- 
mations. Correlation effects are found to be of comparable magnitude as the direct exchange. We 
give quantitative results for realistic quantum dot geometries and find the largest couplings for one 
dimensional systems. 
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Control of spins in semiconductors has been intensively 
investigated in recent years due to its potential for appli- 
cations in spintronics and quantum computation^. Co- 
herent coupling between localized spins is particularly 
sought after because it is a key requirement in any pro- 
posal for spin-based implementation of quantum com- 
putation. Several coupling mechanisms have been pro- 
posed to construct quantum gates between spins in quan- 
tum dots (QDs). These include direct wavefunction over- 
lap using electric gates where the interdot separation is 
smail^ and exchange of a cavity photon mode between 
spins in QDs for a large interdot separation 3 . 

Recently Piermarocchi et at proposed to use an indi- 
rect mechanism to couple the two QD spins at interme- 
diate interdot separations^. Here the interaction is medi- 
ated by virtual delocalized carrier excitations in the host 
material. The virtual excitations are driven by an in- 
terband off-resonance laser that provides optical control 
over the interaction and serves to reduce the bandgap en- 
ergy, thus extending the interaction range. The proposed 
scheme has the advantages of ultrafast optical control and 
long spin coherence times due to the virtual nature of the 
excitations. Combined with the proposal to use Raman 
optical transitions via intermediate trion states for single 
qubit operations 5 , this mechanism provides the necessary 
set of universal gates for quantum computing. This op- 
tically induced indirect spin exchange is a variant of sev- 
eral analogous mechanisms. These include the RKKY in- 
teraction in metals^ and Bloembergen-Rowland coupling 
in direct-gap semiconductors^, the superexchange medi- 
ated by two holes in diluted magnetic semiconductors 
(DMS)&, the magnetic exchange mediated by bound cor- 
related states (excitons)2ii£, and ferromagnetism induced 
by virtual Mn acceptor level-valence band transitions in 
DMS materials^. 

A key ingredient in all of these indirect spin coupling 
mechanisms is the exchange interaction of a localized spin 
with the mediating itinerant excitation. The purpose of 
the present work is to introduce a microscopic formula- 
tion that provides a quantitative description of this ex- 
change interaction by taking into account the effects of 
hybridization of continuum and dot states and the double 
occupancy in the dot. 



For the case of the optically induced indirect interac- 
tion between spins in quantum dots, the spin-spin cou- 
pling is obtained by considering the self-energy correction 
in the continuum electron propagator due to its Coulomb 
interaction with each of the localized spins within second 
order perturbation theory*. The result is a Heisenberg 
Hamiltonian involving the localized spins, with an effec- 
tive positive exchange constant that is given by^ 
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where R is the distance between the dot centers, 8 is the 
detuning of the laser with respect to the electron-hole 
continuum, f2 is the Rabi energy for the light coupling to 
the electron-hole pair and \i is its reduced mass. </(k, k') 
is the exchange interaction between the electron spin in 
the quantum dot and the itinerant electron. 

To calculate J(k, k') we consider a Hamiltonian that 
includes the kinetic energy, the dot potential relative to 
the host material and the electron-electron Coulomb in- 
teractions: 



TL = Ti-o + TLm + Wi 
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Here c\ a (cj^.) is the creation operator for a localized 
(itinerant) electron, n a — c dt7 Cda, and the last term in 
Eq. I|3af) is the on-site Coulomb repulsion. TLm rep- 
resents the hybridization of the localized and itinerant 
electrons, where we include a population dependent hy- 
bridization (2nd term in Eq. I|3b|) h which was absent 
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in previous works on coupling of localized and itiner- 
ant spins. We show below that this latter term makes 
an important contribution to the spin exchange inter- 
action. TL\ contains the spin-independent and spin ex- 
change Coulomb scattering; the latter gives rise to the 
Heitler-London exchange contribution^. The last term 
in Eq. (|5c|l describes the effect of localized and contin- 
uum state mixing. We have neglected scattering pro- 
cesses between carriers in the continuum since the corre- 
sponding effects are not relevant to the problem we wish 
to solve. V k = J drVdotVkl 1 ")^! 1 ") i s the tunnelling am- 
plitude, where v?k(?~) [VdW] is the itinerant (localized) 
electron wave function, also used to calculate he various 
Coulomb integrals in Eq. (J2J. 

We aim at bringing the Hamiltonian TU = TLo + TLm 
to a form similar to that of an s-d spin exchange Hamil- 
tonian by applying a canonical transformation 
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The unitary operator S is constructed to eliminate TLm 
to first order by requiring TLm = — [S, TLo] and is given 
by 



kcr 



h.c. 



where 



ctk = 



V k +T k 



U + E d 



E k 



(3 k = 



Ed — Ei. 



(5) 



(6) 



This is a generalized form of the Schrieffer- Wolff transfor- 
mation, which was first used to establish the connection 
between the Anderson and Kondo models 12 . It produces 
a contribution to the spin exchange arising from the cor- 
relation and hybridization terms in TL', which is given to 
first order by 

jW(k,k') = ~[f3 k V k *, -a k {V k , + T k ,)*} + [k~k'}* . (7) 

This contribution vanishes when correlation effects are 
neglected (U, T k — > 0). We find that the first order 
result, given in Eq. (Q, is inadequate because it requires 
that TLm would be a small perturbation to TLo, which is 
not the case generally. It is therefore necessary to sum 
up the infinite series in the transformed Hamiltonian 
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where [S,H M ]n = [S, [S, [S,TL M ]--}}- To this end we 
use a method suggested by Chan and Gulacsii^ but em- 
ploy a different strategy to solve the problem. The first 
term in the series in Eq. JBJ is 
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where J\{k,k') — 2J^\k,k') is given in Eq. JTJl, 
and the rest of the coefficients in Eq. @ are 
Pi (*;,*;') = \[a k V£ -!3k{V k - +T k ,)*} + [k <-> k']*, 
#i(M') = (3 k V k * + (3* k ,V k , Gi = 2E*PkV k *, and h = 
2 J2k [otk{V k + T k )* - hkV k \- Calculating the second or- 
der term we find that it has the same form of TLm apart 
from higher order correlation terms. We estimate the 
magnitude of these continuum scattering terms by ne- 
glecting off-diagonal contributions and placing lower and 
upper bounds on the occupation numbers. This proce- 
dure brings all higher odd orders in the series to the form 
of Eq. @ , and we are able to sum the series by solving 
the following set of recursion relations for the several co- 
efficients 
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Equations (|1U[) are obtained from the lower bound in the 
higher order contributions, and a second set of equations 
is obtained for the upper bound case. 

The exchange contribution is obtained from the odd 
orders of the series, which also contain additional terms 
that renormalize the original Hamiltonian l|8a[l. The 
even orders also are summed up and renormalize the hy- 
bridization Hamiltonian l|3b(l . Figure la shows the result 
of the series summation for the diagonal part of the ex- 
change, J(k, k). Since it differs appreciably from the first 
order result of Eq. (JJJ) , the residual hybridization in the 
even orders need not be small, as seen in figure lb. Thus 
we need to perform a second canonical transformation 
that is defined by applying Eqs. (|5I6|I to our renormal- 
ized Hamiltonian. This second transformation eliminates 
the next order in the hybridization terms and further cor- 
rects the resulting exchange contribution. The process is 
reiterated until we fully eliminate the hybridization part 
of the Hamiltonian, as shown in figure lb. 

This procedure of applying a set of nested Schrieffer- 
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FIG. 1: (color online) (a) Diagonal matrix elements of the 
spin exchange interaction between localized and itinerant elec- 
trons in a 2D host. The figure presents the first order result 
from Eq. Q (dashed line) , intermediate results after the sum- 
mation in Eq. (JgJ (dotted lines), and final results obtained 
by performing a set of transformations (solid lines); (b) The 
corresponding hybridization terms: tunnelling amplitude, Vk, 
and population-dependent hybridization term, Tj,. 



Wolff transformations is essential to obtain quantitative 
results for the kinetic exchange interaction, which is the 
one that results from the hybridization terms in Eq. l|3bjl . 
As seen in Figure la, J(k, k') is ferromagnetic after one 
transformation, which would differ from other results for 
this kinetic exchange contribution, e.g., in a renormaliza- 
tion group approach^. Only after a set of transforma- 
tions (typically 10-20) are the renormalized hybridization 
terms eliminated and the antifcrromagnetic nature of the 
interaction is restored, albeit with a modified magnitude 
compared to the first order result. The results calculated 
with the lower and upper bounds discussed above are re- 
markably close to one another. We have verified that they 
coincide within 10% for a wide range of geometries and 
dot potentials. In order to address the possibility that 
the off-diagonal contributions from the continuum scat- 
tering terms might alter our results, we have estimated 
them from limits of the off diagonal density factors. The 
results lie between the lower and upper bounds given in 
figure la; thus we believe that our summation represents 
the complete Schrieffer- Wolff transformation with a good 
accuracy. 

Since the kinetic exchange interaction that we calculate 
from the transformed Ti! is antifcrromagnetic, it com- 
petes with the ferromagnetic exchange given by the sec- 
ond term in Eq. (|3cfl . Thus, an accurate evaluation of the 
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FIG. 2: (color online) (a) Optically induced spin-spin ex- 
change coupling in a 2D host vs. the distance between the 
centers of the dots for dot radius Rd = lOnm, potential height 
V e = 150meV and several values of laser detunings. Right 
axis shows the coupling values including excitonic effects; (b) 
Same as a for a quasi ID host and cylindrical dots with 
L z — lOnm, Rd = 5nm and potential height V e = 80meV. 



former is important as it can lead to an order of magni- 
tude difference or even a change of sign in the total spin 
exchange coupling between localized and itinerant elec- 
trons. A full transformation is also valuable in the case 
where U+Ed > 0, leading to a divergence of in Eq. JJjJ. 
Here, the kinetic exchange contribution is dominant and 
cannot be obtained via a perturbative approach. This 
regime corresponds to dots with small size (Rd < 5nm) 
and shallow potential (barrier < 80meV), which are not 
typical for physical systems and are not considered here. 

In figure 2 we show the results for the spin-spin cou- 
pling J12 [Eq. Q], incorporating all the exchange con- 
tributions. Figure 2a shows the spin coupling for lateral 
cylindrical dots in a two-dimensional quantum well. The 
results for vertically stacked cylindrical dots in a quasi 
one-dimensional wire are given in figure 2b. Here we 
used m e = 0.07m, mh = 0.5m and f2 = O.lmeV. The 
localized electron wave function was taken to be a com- 
bination of Bessel functions in the lateral direction and a 
combination of Cosine and Exponential functions in the 
z direction. It is seen that the spin coupling is more than 
an order of magnitude larger for the one-dimensional case 
than for the two-dimensional case. 

The Coulomb interaction between the intermediate vir- 
tual electrons and holes results in an enhancement of the 
oscillator strength at the optical and spin vertices due 
to the exciton wave functions^. We have evaluated the 
dominant contribution of the electron-hole interaction to 
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FIG. 3: (color online) (a) Spin-spin coupling in a 2D host at 
a dot separation of 21nm vs. the dot potential for two dot 
radii. The laser detuning is 8 — 0.5meV. Right axis shows 
the coupling values after excitonic corrections; (b) Same as a 
for a quasi ID host and cylindrical dots with Rd = 5nm and 
several dot heights. 



J\2- It results in an increase of up to two orders of mag- 
nitude in the two-dimensional case and roughly one order 
of magnitude in the one-dimensional case (right axes in 
figure 2). Thus, the excitonic effects reduce the difference 



in J i2 between the two geometries. 

Figure 3 shows the dependence of the spin-spin cou- 
pling on the dot potential and size. Larger dots give 
larger couplings but necessitate larger separation to avoid 
overlap. The increase in the coupling as the dot poten- 
tial decreases is mainly due to the reduction of the kinetic 
exchange contribution. 

A technologically viable way to increase J12 is by using 
a microcavity. This can be done by growing distributed 
Bragg Reflector layers on the top and bottom of the ac- 
tive semiconductor layer containing the QDs. Placing the 
active layer at the antinode of the microcavity increases 
the electric field by orders of magnitude, and thus in- 
creases the Rabi energy at the optical vertices in Eq. . 

We have shown that the effect of hybridization of con- 
tinuum and dot states produces a sizable contribution 
to the exchange coupling between localized and itiner- 
ant electrons. For certain dot geometries this kinetic 
exchange can even lead to a change of sign in the spin 
exchange interaction. A set of canonical transformations 
with summations over higher order terms provides a use- 
ful tool to evaluate the spin exchange interaction. Our 
transformation of the Hamiltonian (0) captures the mul- 
tiple scattering processes involved in the interaction be- 
tween the localized and itinerant carriers, and it provides 
the first microscopic description that accounts quantita- 
tively for the exchange interaction 16 . Our formulation is 
also applicable to other systems of localized spins cou- 
pled by carriers, such as electrons bound to donors^, 
magnetic impurities^ and nuclear spins^. 
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